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Abstract: We present a general analysis of B s d — > £ + £~ and Kl — > £ + £~ decays 
in supersymmetric models with non-minimal sources of flavour mixing. In spite 
of the existing constraints on off-diagonal squark mass terms, these modes could 
still receive sizeable corrections, mainly because of Higgs-mediated FCNCs arising 
at large tan (3. The severe limits on scenarios with large tan (3 and non-negligible 
d l R{L,)~d 3 R(L) mixing imposed by the present experimental bounds on these modes and 
AB = 2 observables are discussed in detail. In particular, we show that scalar- 
current contributions to Kl — > £ + £~ and B-B mixing set non-trivial constraints on 
the possibility that B s — > £ + £~ and Bj — > £ + £~ receive large corrections. 
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1. Introduction 

The rare dilepton decays B s ^ — > and A^ — ► ^ + ^~ are ideal candidates to deeply 
investigate the dynamics of quark-flavour mixing. On the one side, they are highly 
sensitive to physics beyond the Standard Model (SM), like all processes mediated by 
flavour- changing neutral-current (FCNC) amplitudes. On the other side, the helicity 
suppression of vector contributions gives to these modes a unique sensitivity to scalar 
currents, which are typically negligible in other rare decays. 

In this paper we analyse the rare dilepton decays in supersymmetric (SUSY) 
models with flavour non-universal soft-breaking terms. As is well known, this class 
of models is strongly constrained by currently available data on AF = 2 and AF = 1 
amplitudes (see e.g. Ref. [0]). However, a large fraction of the parameter space still 
remains to be explored. The only conclusion that can be drawn at present is that 
the possible new sources of flavour mixing, if any, must have a rather non-trivial 
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structure. For instance, as emphasized recently in Ref. @], the possibility of a large 
bR-SR mixing, motivated by SO(10) unification and neutrino data, is still open. 

Supersymmetric contributions to rare FCNC decays have been widely discussed 
in the literature. However, most of the analyses are focused on scenarios with tan f3 = 
0(1). As we shall show, if tan/3 = 0(1) there is little to learn from the dilepton 
modes, which cannot be extracted elsewhere. On the other hand, things change 
substantially at large tan/3, because of the scalar FCNCs induced by the effective 
coupling of to down- type quarks (the so-called non-holomorphic terms) || Q . 

As discussed in Refs. [[|-@, at large tan/3 the rates of B s ^ — > £ + £~ decays 
could be substantially enhanced with respect to the SM even without new sources of 
flavour mixing. As pointed out in 01] , the same conclusion does not hold for — > 
£ + £~. In this paper we analyse how the situation is modified by non-minimal flavour 
structures. Interestingly enough, in this case also Ki — > £ + £~ play an important 
role, setting stringent bounds on off-diagonal squark mass terms. 

The issue of SUSY contributions to B s d — > £ + £~ decays in the presence of non- 
universal squark mass terms has already been addressed in Refs. J6], [/J. Here we 
extend these works by discussing all possible tan /3-enhanced terms, both those arising 
at higher orders and those not related to scalar FCNCs. We also clarify the issue 
of the diagonalization of the quark mass matrix in the presence of non-universal 
non-holomorphic terms. Last but not least, in our analysis we take into account the 
constraints from scalar-current contributions to Kl — > £ + £~ and B-B mixing, which 
have not been discussed before. 



2. Generalities 

The part of the AF = 1 effective Hamiltonian relevant to P\d l dP\ — ► decays is 
C 2 M 2 

H eft = F_W [Ca0a + C < ry + CsQs + CpQp + Cs&s + C , p(D , p] + kc-) (2.1) 

where 

O s = d R di££ , 0' s = d\d? R ££ , 

P = d R d{I l5 £ , 0' P = d l L d R £ lb £ . (2.2) 

According to this normalization of the effective operators, writing hadronic matrix 
elements of axial and pseudoscalar currents as 

M 2 

mi,l^\P{p)) = ipjp , (0\qh 5 q j \P(p)) = -ifp 7 / , , (2-3) 
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the most general expression of the short- distance (SD) contribution to P[d % d^\ —>■ 
branching ratios reads 



B SD (P 



^4 71^4 

-JpMpT P m e 



Amj 

Ttp~ 



1/2 



X 



Amj 
~Mp~ 



M P (C S 



C' s ) 



2me(mi + rrij) 



M 2 P (Cp - C'p 



2me(rrii + rrij] 



+ Ca — C' A 



(2.4) 



where the Wilson coefficients are understood to be evaluated at a scale of O(Mp). 
In B s d modes this is the only relevant contribution to the branching ratios, whereas 
in K L decays also long-distance effects have to be taken into account ||. 
Within the SM only Oa has a non- vanishing Wilson coefficient, 



cT(m 2 w ) 



m t 



(2.5) 



where the charm contribution, which receives sizeable corrections in the running 
down to the meson scale, plays a non- negligible role only in the Kl case . In the 
SUSY scenarios we are considering in this paper, all the Cj could have, in principle, 
a relevant impact on the decay amplitudes. 

We stress that, since the ratios {Cs,p — C' SP )/me are independent of me, to a first 
approximation the relative weight of the various contributions is independent of me, 
whereas the overall branching ratios scale like m\. This implies that i) e + e~ modes 
are never interesting; ii) the difficult detection of B(B s d — ► r + r~) with respect to 
B(B S) d — > /U + /i~) is partially compensated by an enhancement factor ~ (m T /m M ) 2 . 



3. Axial operators 

The contributions to the axial operators Oa and 0' A can be divided into two big 
categories: Z penguins and box diagrams. Those in the former class receive con- 
tributions from charginos, neutralinos and gluinos, while gluino exchange does not 
contribute to the latter class. 

As discussed by several authors [EI]— [13], box diagrams always play a subdom- 



inant role: in the case of q^-Qji mixing the Z penguins are by far dominant, in the 
case of S'?7(2)x / -conserving mixing {q\-q 3 L or q l R — cp R ) the constraints from AF = 2 
processes forbid sizeable effects. 

Z-penguin diagrams correspond to dimension-six effective operators of the type 
(f^P^q^PjH^ D , where H denotes the Higgs field. For this reason, to be non- 
vanishing, these diagrams require a double chirality flip, which can arise either by 
the q % L-q R mixing or by the gaugino-Higgsino mixing (for charginos and neutralinos 
only). Also in this case, AF = 2 constraints forbid sizeable effects in the case of 
S , ?7(2)x-conserving squark mixing. As a consequence, possible large effects (with 
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respect to the SM) are controlled only by the q l L-<i R mixing terms, both in the up- 
and in the down-type sector. 

3.1 Charginos 

Involving up-type squarks, chargino contributions to Z penguins are mildly sensitive 
to the value of tan j3. These contributions, which we report here for completeness, 
have been widely discussed in the literature in the context of B(K) — > X s (X d )£ + £~ 
processes |II]]-Jl4|]. Before any mass expansion, we can write [ [3~2|j 



\- Ca ^X ~ ^ L il L mk F nmlk (3.1) 



where 



U mk = -g{t u )\ Lk V mX + y-V^Y% Rk V m2 , 

''nmlk = ^nl^lm k(x m k, X nk ) 2C/ m lC/^ n 5[k y/ % mk% nkj (% mk-i -^nk) 



-5 mn fl(T u )l Tl k(x mk ,x lk ) ; (3.2) 



g is the SU(2) i gauge coupling, V denotes the CKM matrix, = \r2m Uq /(v sin (3) = 
gmu q / (V%M W sin/?) and the loop functions k(x,y) and f(x,y) can be found in the 
appendix. 1 Here V and U are the unitary matrices which diagonalize the chargino 
mass matrix [U*M X V^ = diag(M Xl , M X2 )] and T u is the one which diagonalizes the 
up-squark mass matrix, written in the basis where the d % L — u° L — \n coupling is family 
diagonal and the d\ — u J R — Xn one is ruled by the CKM matrix. 

The contribution of chargino Z-penguin diagrams to C' A necessarily involves two 
down-type Yukawa couplings (yfyj), which cannot be both of the third family. As a 
result, even at large tan/?, this contribution is negligible with respect to the one in 
Eq. (|3). 

A detailed discussion on the relative weight of the various terms appearing in 



Eq. (|3.1| ), up to the second order in the mass expansion, can be found in JT2|, |14| . Due 
to the off-diagonal CKM structure, the zeroth order term is not vanishing. However, 
this is not very interesting since it reaches at most 0(10%) of the SM contribution: 
the potentially largest effects are induced by terms of first or second order in the 



expansion of the squark mass matrix around the diagonal [0, A substantial 
simplification arises if, in addition to the squark mass matrix, also the chargino mass 
matrix is perturbatively expanded around the diagonal. In this case we can write 



[C 



Alv° 



(jSM 



V 3j V 3i V 3i V 3j 



+ ..., 

(3-3) 

1 As usual x a p — M^/Mg, with the indices (n, m) and (/, k) denoting chargino and squark-mass 
eigenstates, respectively. 
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where 

Vrl)v ~ (M | } > d - , (3.4) 

and the dots denote terms not enhanced by inverse CKM matrix elements, which 
can be safely neglected. The explicit expressions of the adimensional functions /i 2, 
which depend on the various sparticle mass ratios (and mildly on tan/5), can be 
extracted from the full result in Eq. ( PI) : for M~ tR /Mu L > 1/2 and M Xn /M Uh > 1/3, 
one finds |/i| < 0.1 and I/2I ~ 0.4 (the upper figures are obtained for the minimal 
value of Mi ). 

As anticipated, large effects are in principle possible in the presence of large 
u l L -u J R mixing. However, both in the kaon and in the B system this possibility 
has recently been severely restricted by new experimental data on the non-helicity- 
suppressed modes. In the kaon case a clean and stringent bound is set by B(K + — »■ 
n + uu) Similarly, in the B system severe bounds are obtained by the recent data 



on both exclusive and inclusive Bd — > X s £ + £ transitions [fL6 j. As a consequence of 
these new data, we are led to the model-independent conclusion that \Ca/Ca\ < 2 



in the s — > d channel |17| and is even closer to unity in the b — ► s case ||18|| . The 
only channel without strong experimental bounds is the b — d one. However, once 
the vacuum-stability bounds of Ref. Jl9| are imposed on ( 5^)31, also in this case one 
can at most obtain \Ca/C a m \ = 0(1). 

3.2 Gluinos 

Contributions to Z penguins which involve down-type squarks, with gluinos and 
neutralinos, are very sensitive to the value of tan (3. In particular, at small tan (3 
their effect is completely negligible because of the smallness of left-right down-type 
mass insertions. This conclusion is not necessarily true at large tan/3, where, in 
principle, left -right down-type mass insertions can become as large as up-type ones. 

The largest effects are expected from gluino exchange, whose contribution to 
before any mass expansion, can be written as 

[CAT -^rAAkhalk, (3.5) 



89 



where 



4 = -V2g s T b tg , 
A\ k = -V2g s T a (T D )l Lk 



haik = S ab T kQR (T y qRl k(x gh x kl ) , (3.6) 

g s and T a denote SU(3) C coupling and generators, and T D is the matrix that diag- 
onalizes the down-squark mass matrix (in the super-CKM basis). Contrary to the 
chargino exchange, in this case there is no zeroth order term in the mass expansion 
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and C' A is not negligible a priori: [C' A \g = [Ca]^ (ql Qr)- However, the bounds 
on down-type mass insertions obtained from AF = 2 processes are much stronger 
than in the up-type sector (because of the presence of AF = 2 gluino-box ampli- 
tudes) U). Once these constraints (and the vacuum-stability ones) are taken into 
account, it is easy to realize that in all types of d l — > d? transitions both [C' A \g° and 
\C A \ Z g can barely reach the magnitude of |C| M |. This conclusion is important since 
it ensures that the model-independent constraints on Ca discussed before cannot be 
invalidated by an accidental cancellation between C' A and Ca- 

The whole discussion on axial-vector operators can thus be summarized as fol- 
lows: (i) at present, in all types of d % — > d? transitions, axial- vector contributions 
can reach the SM level at most; (ii) the most efficient way to constrain (or detect) 
non-standard axial-vector contributions is by means of non-helicity-suppressed decay 
modes. 



4. Scalar and pseudoscalar operators 

As discussed in Refs. @-||, in the limit of large tan (3 Higgs-mediated scalar FCNCs 
can strongly modify, or even overcome, the SM amplitudes of rare helicity-suppressed 
modes. The origin of this effect is related to the appearance, at the one-loop level, 
of an effective coupling between Hu and down-type quarks f2"0fl . As discussed in J3J 



see also ||22||), a necessary condition to take into account all tan /3-enhanced terms in 



FCNC amplitudes is the diagonalization of the down-type quark mass matrix in the 
presence of these effective couplings. Here we generalize this approach to the case 
where the effective d l R d 3 L H^j coupling contains also non-minimal sources of flavour 
mixing. 2 

Restricting the attention to the down-type neutral-Higgs sector, the effective 
down-type Yukawa interaction we need to consider is given by 

£f = d R [yfSijK + EuH*] d{ + h.c. , (4.1) 

where Eij parametrizes the coupling generated at one loop by both chargino and 
gluino diagrams. Here we do not make any specific assumption about the flavour 
structure of E^, but require that the off-diagonal entries can be treated as pertur- 



bations with respect to the (full) diagonal mass terms in fl4.1|) . Then proceeding as 



2 We stress that the so-called resummation of the tan /3-enhanced terms, in the limit where 
corrections of 0(M^/M 2 ) are neglected, is nothing but the diagonalization of the effective Yukawa 
interaction, which includes the dimension-four effective operators appearing only at the quantum 
level |Q, [22) . In particular, the computation of the coefficients of all relevant dimension- four effective 
operators to order (ai/ir) 1 , together with the diagonalization of the effective Yukawa interaction, 
leads to a resummation of all terms of order (aj tan f3/ir) n in the limit M{y <C M 2 . 

3 In particular, this means we ignore the fine-tuned scenario where En tan/3 w —yf. 
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in Ref. [Q with a perturbative diagonalization of the mass terms, we find 



eff {i>j) _ yf tan g 
^-fcnc - -d 



y[_y[ 

r.d „.d 



1 ;H°;-H° d 



tan/? M 



« + t-c, (4.2) 



where 

yf = yt + E ii tasil3 . (4.3) 



The analogue of (|4.2j) for z < j is obtained with the exchange i <-> j in all indices of 
Yukawa couplings (both and |/ d ), but not in the £7^ terms. 

In the absence of non-minimal sources of flavour mixing, one has 

En = El + El , E% = e g yf6 lv E* = eyyfy^V", (4.4) 

where V° denotes the bare CKM matrix (without the ey tan f3 corrections induced 
by re-diagonalization) and 

e 9 = f ( X gdL -i x d H d L ) 5 e Y — ^ Q~2 f ( X /J-u L j x u R UL ) • (4-5) 

d L u L 



Using Eqs. ( |4.2| )-(fOD, it is easy to recover the results of Refs. H |J. We shall now 

% and £^ 



proceed with the evaluation of Ef- and Ef- in the presence of non-minimal sources of 



flavour mixing. 
4.1 Charginos 

The Eij term induced by chargino-up-squark loops can be written as 

E% = ^- 2 V nl KkG n ik , (4.6) 



where 



Kk = yfi^tMn, (4.7) 

G nlh = y t rl L (AuV )l tR (r u )l Rl ^f(x lk , x nk ) . (4.8) 



As usual, the trilinear soft-breaking term has been decomposed as Ajj x Yjj, and we 



have neglected all up-type Yukawa couplings but the top one. Since oc yf, only 



the case i > j is relevant, and it is also clear that the E*^ term in Eq. (|4.2j ) can safely 
be neglected. 

Integrating out the heavy neutral Higgs fields, in the large-tan/? limit, we then 
obtain 

\ r ]h°,H°,A° = \n ]h°,H°,A° 71-2 yjEfjt&u-P 

[ Slx " 1 Plx " IGlM^Ml [1 + tan/3(e 9 + e Y y?5 i3 )} ' ^ 



-7- 



while the coefficients [C^p]^ '^ '^, suppressed by the factor y^/yf, turns out to be 
negligible. Due to the off-diagonal entries of the CKM matrix, even with a diagonal 
up-squark mass matrix Ef,j 7^ 0. This zeroth order term leads to |4|: 



A*, tan 3 B n.Af (v.... ,■ .. \ 

(4.10) 



m dt m e m 2 A^-tan 3 /3 /j,Af (x^, x UrUl ] 



4M^ [1 + tan/3(e s + e Y y 2 5 i3 )} M^M^l + e 9 tan/3] 

where, expressing V° in terms of the physical CKM matrix, 

Vi 3 V 3j (i = 3) 

"l+l 

l+e g tan /3 



X i>3 



v 3 *v 3j 



Interestingly, this zeroth order term in the mass- insertion approximation (MIA), 
turns out to be the most interesting contribution induced by chargino exchange. 
Indeed, the only possibility to remove the CKM suppression factor V 3 j (the only one 
relevant in B decays) is by means of a double mass insertion: one in the chargino 
line, the other in the Z-squark line. This leads to a potential enhancement factor 
(8rl)3j8 X /^3i> which is partially compensated by the suppression of the loop function 



with two mass insertions, i.e. a structure very similar to Eq. (3J3). Then, taking 
into account the bounds on (5^ L ) 3 j discussed in Section pTT] , it turns out that this 
contribution can be at most of the same order as the zeroth-order term, or the 
minimal-flavour-violating (MFV) term. 

The suppression factor V 3 *, which is relevant only to the kaon case, can be more 
easily removed: we simply need an off-diagonal trilinear term not proportional to V^*. 
This leads to a net enhancement factor {S^jijis/V^, which could reach one order of 
magnitude. However, as discussed in Ref. Q, the zeroth-order chargino contribution 
to the scalar — > amplitude amounts only to a few per cent of the SM one. 
Thus this potential enhancement is not particularly relevant. 

4.2 Gluinos 

The Eij term induced by gluino-down-squark loops can be written as 

% = ^AlKkhalk , (4-12) 



where 



Bl k = Ai k (i L ^ i R ) , (4.13) 



di 



J ba ik = ^ty q {V D )\ Lll ^f{x gdv x dkdl ) . (4.14) 



As can be noted, the gluino contribution to E^j has a structure rather different from 
the chargino one. The latter is different from zero already at the lowest order in the 
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MIA, but it is necessarily proportional to the Yukawa coupling yf. On the other 
hand, Ef,, vanishes at the zeroth order in the MIA; however, with a suitable mass 
insertion it can be made proportional to any of the down-type Yukawa couplings. As 
a result, in the gluino case also the C' s P coefficients can play an important role. 

Assuming for simplicity i > j, the leading contributions to the scalar operators 
induced by Ef^ are given by 

[Sla ~ [Pla ~3g 2 M%>m [l + tan/?(e ff + e y ^ 3 )] 2 



x 



(^LL,)ij X dd L f( X 9d L j x d R d L ) 3j x dd L f( X gd L > X d R d L , Xd R d L , 1) 



(4.15) 

< r tih°H°A° = \ r i ih°H°A° = ^9s_ m di minMg tan 3 (3 

1 - ' 1 H " W M \M [l + tan/3(e g + e y y t 2 5 i3 )] 2 



d 

D 2 Ufa D D 

ij X dd R fi X gd R i x d L d R i 1) + d^LL)i^RR)^3 X dd R J "i x gd R i x d L d R , %d L d R 

Vi 

(4.16) 



where the off-diagonal elements in the squark mass matrix have been parametrized 
in terms of flavour-changing insertions in the sfermion propagators, defined as in 
Eq. d3.4|). The terms arising at the second order in the MIA are relevant, and play a 



very important role, only in the s — > d case [(Slurr))^®]- in deed, similarly to what 
happens in chargino- induced Z penguins 0], only with a double flavour-mixing of 
the type (2 — ► 3) x (3 — ► 1) can the large Yukawa coupling of the third generation 
contribute to s — > d amplitudes. 

Before analysing the phenomenological implications of these formulae, we point 
out some differences with respect to previous analyses: 



Contrary to Refs. 0, the expressions ( j4.15| ) and ( [4.1 6|) contain the leading 



(flavour-diagonal) etan (3 corrections to all orders. As pointed out in [|]], this is 
a numerically important contribution for tan/3 ~ 50, where the largest effects 
are expected. In the limit where higher-order etan/3 corrections are neglected, 
we find that C 9 S p is proportional to tan 3 /3, similarly to Cg p . This behaviour, 
which is explicit in Ref. 0, is somehow hidden in Eqs. (4.4)-(4.5) of Ref. 
(where the effective Yukawa interaction is not diagonalized) . However, the 
tan 3 f3 behaviour become manifest also in Ref. J7| after an expansion of the T D 
matrices: in this limit we fully agree with the overall normalization of Ref. [|7|]. 4 

Within the limiting factors associated to the MIA, our formulae contain as a 
special case the scenario C of Ref. M. In our opinion, such a scenario should 



4 This normalization differs from the one of Ref. ||, which contains an extra factor of 4. Wc 
thank J. Urban for a clarifying discussion about this point. 
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not be denoted clS db CclSG of minimal flavour violation. Indeed diagonal - 
but flavour-non-universal — squark mass matrices do provide a new source of 
breaking of the SU(3) 5 flavour group |22|] and indeed they break the SM 



relation between Yukawa matrices and FCNC amplitudes. Since the purpose 
of this analysis is to identify the potentially largest effects generated by non- 
minimal sources of flavor mixing, and to identify the main differences with 
respect to MFV scenario (defined as in Ref ]22|]), we prefer to treat all non- 
minimal cases in the same way, using the MIA. For the same reason, we shall 
not consider subleading terms such as the neutralino amplitude, which become 
important only in presence of cancellations among the leading terms M. 



5. Phenomenological bounds 

As discussed in the previous two sections, only scalar (and pseudoscalar) operators 
can induce large enhancements of B d , s — > £ + £~ rates over the SM expectations. 
Among the contributions to these operators, the one induced by chargino exchange 
is likely to be dominated by the MFV term, widely discussed in the literature. On the 
other hand, the gluino contribution is extremely sensitive to possible new sources of 
flavour mixing. For this reason, in the following we shall concentrate only on gluino- 
mediated scalar amplitudes. 

5.1 B diS -> 

In the limit where tan/3 is large (tan/3 ^ 30), Ma is rather light (Ma ~ 500 GeV) 
and the off-diagonal mixing terms (Sll<rr))^ are no ^ vanishing, the gluino (scalar) 
amplitude could dominate both B d — > £ + £~ and B s — > £ + £~ rates. To get a rough 
idea of the possible effects, in the limit of degenerate sparticles we can write: 



-< B -^> ^,^(^Y 1( S UrMt ML,- (5.2) 



B SM (B S fi+fi-) \M A J g + 

These approximate numerical expressions hold at first order in the MIA, in the limit 
where we neglect any interference between supersymmetric and SM amplitudes and 
between supersymmetric contributions with different mass insertions. Moreover, we 
have neglected ey^f with respect to e g in the denominators of Eqs. ( |4.15[ ) and ( |4.16| ), 
and we have set e g to the value it assumes in the limit of degenerate supersymmetric 
scales, for /i > (e g ~ 0.012). 5 A more complete illustration of the full dependence 
of B(B d — > fi + fi~) from various supersymmetric parameters is provided by Fig. ||. 

5 The sign of [i is chosen according to the indications from b — > sj and (g — 2)„ [p3[. 
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Figure 1: Left: B{Bd — > /U + fi ) normalized to the SM prediction as a, function of x ^cLl •> 
for (<5f L ) 3 i = 0.1, M A = 200 GeV and different values of tan/3 and x gd L - Right: present 
bounds on \{^ii^ RR ))?,i\ imposed by the (scalar) gluino contribution to B{Bd — > /U + ^x~) as 
a function of tan/? and M^, for x^l = x gd L = 1- 



Using the present experimental bounds 

^(Srf -> < 2.0 x 10~ 7 90% C.L. m, (5.3) 



B (S, -> /x + //~) < 2.0 x 10~ 6 90% C.L. |25j , (5.4) 



we can derive a series of constraints on the squark flavour-mixing terms appearing 
in Eqs. Q4.15| ) and ( |4.16| ). These are reported in Table |I| and illustrated in Fig. |1] (we 
show explicitly only the 1 <-> 3 case since the 2 *-> 3 is completely analogous, except 
for the rescaling of the SM contribution). All the bounds have been obtained in the 
limit of complete degeneracy among the eigenvalues of down-type squark masses. 
Moreover, we stress once more that these bounds are valid only in the limit where 
the gluino scalar amplitude corresponding to a given flavour-mixing term dominates 
completely the decay rate. If (8^)31 and (5 RR )3i are of the same order of magnitude, 
interference effects can no longer be neglected and leads to 0(1) modifications of 
these bounds. For instance for (S^ L ) 3i = (S RR ) 3i the bounds become more stringent 
by about a factor a/2. 

As already pointed out in Ref. |5[], for sufficiently large values of tan (5 the bounds 
on \(&ll(rr))3i\ are more severe than those derived from the gluino-box contribution 



to Bd-Bd mixing |26|]. However, we find that the effect is substantially reduced, with 
respect to what discussed in Ref. |J, by the higher-order etan/3 terms. 

Already for tan/3 = 30, the upper bounds obtained on \(8 RR ( LL j) 32 \ are much 
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\($LL(RR)hl\ 


\( S LL(RR))&\ 




x gd L 


tan /3 = 30 


tan/3 = 50 


tan (3 = 30 


tan (3 = 50 


1 


1 


0.089 


0.026 


0.22 


0.065 


0.25 


1 


0.13 


0.035 


0.33 


0.086 


0.25 


4 


0.16 


0.043 


0.39 


0.11 




\ Re ($LL(RR))^\ 




)23(5 LL(jRR) ) 3 l]| 


x /id L 


x gd L 


tan (3 = 30 


tan/3 = 50 


tan /3 = 30 


tan /3 = 50 


1 


1 


0.017 


5.0xl0~ 3 


0.7xl0~ 3 


2.0xl0~ 4 


0.25 


1 


0.025 


6.6xl0~ 3 


l.OxlO- 3 


2.7xl0" 4 


0.25 


4 


0.030 


8.0xl0- 3 


1.4x10-3 


3.7xl0" 4 



Table 1: Present constraints on the off-diagonal squark mass terms that rule the (scalar) 
gluino contribution to B s d — > fi + n~ (upper) and Ki — > /J + /U~ (lower). All the bounds 
have been obtained for M A = 200 GeV and scale as (Ma/200 GeV) 2 . 



more severe than those that could be derived, in the near future, from B s -B s mixing 
or the CP-violating asymmetry in Bd — > 4>K S . As a result, if the scenario with a large 
^r-Sr mixing discussed in Ref. is accompanied by a large value of tan/3, which is a 
natural possibility in SO(10) models, the most interesting observables to look at are 
B (B St d — > ji + n~). The clear signature of this scenario would be B(B S — ► yU + /i~) well 
above the SM expectation and, at the same time, a ratio BiBd — ► fi + fi~)/B(B s — > 
/i + yU~) different from (V^/V^l 2 (i.e. the value expected in the absence of new flavour 



structures |22|| ) 



5.2 if L -> 

Within the SM the short-distance contribution to the — > fi + {i~ branching ratio 
is predicted to be B SM (K L -> fi + fi~) S B = (8.7 ± 3.7) x 10~ 10 [HJ. This expectation 



cannot be directly confronted with experimental data because of the long-distance 
amplitude generated by the two-photon intermediate state. Indeed the precise exper- 
imental measurement of B (Kl — > /i + /U~) |57] turns out to be completely saturated 
by the absorptive two-photon contribution, signaling a cancellation between short- 
distance and dispersive long-distance amplitudes. Taking into account the estimates 
of the dispersive long-distance amplitude in Refs. || , we shall assume in the following 
the conservative upper bound 

B (K L -> /i V)sd < 3.0 x 10- 9 , (5.5) 

which should be regarded as a reference figure rather than as a strict bound. 

As can be read from Eqs. ( J4.15|) and Q4.16I ), gluino contributions to the K L — > 



£ + £~ amplitude suffer, at the first order in the MIA, by a double strong suppres- 
sion: the smallness of the strange Yukawa coupling and the indirect constraints on 
(^ll(rr)) 21 BUI- However, the SM amplitude of Kl — > £ + £~ decay is also affected 
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Figure 2: Constraints on |Re[(5^^ Li ^)23(^fx(^m)3i]| se t by - ► H + H ■ Left: the 
bound as a function of x^cL L , for — 200 GeV and different values of tan {3 and Xgd L . 
Right: the bound as a function of both tan/3 and Ma for x^d L = x gd L = !• 



by a strong suppression due to the CKM hierarchy [A(Kl — > £ + £~) SM ~ A 5 ]. As a 
result, for large values of tan/3 this supersymmetric amplitude could even become 
the dominant short-distance contribution. Adopting the same approximations as in 
Eqs. ( |5.1| ) and (p.2|), we can write 

B°W(K L ^ fi+fj,-) _ /200\ 4 

If |Re(5^ L )i2| 3> |Im(5f > L ) 12 |, as happens in many reasonable scenarios, the bound on 
l-^ e (^LL(_R J R))2i| we can extract from this amplitude (see Table P becomes competi- 
tive, for tan f3 ~ 50, with those derived from the gluino-box contribution to K°-K° 
mixing ||26|| . 

As already emphasized, in rare K decays it is often necessary to go up to the 
second order in the MIA in order to take into account all possible large effects |TjJ. 
This occurs also in this case where, according to Eqs. (|4.15|) and ( |4.16|) , at the second 
order in the MIA we can get rid of both suppression factors affecting the term with 
a single mass insertion. In the usual limit of degenerate supersymmetric scales, we 
can write 

B9W(K L ^ n+fj,-) _ /200\ 4 
B^(K L ^fi+fi-) SB ~ \M A ) 

(5.7) 



Re(5 D 



LL{RR))1l 



3 X 10 



-3 



/ tan p \ b 
I 50 ) 



(5.6) 



^ e [(^RR(LL))^LL(RR))3l\ 



1 X 10" 4 



/ tan g y 
\ 50 ) 



]+!(¥)] 



i4 - 





(*&)sil 


^1(032(^)321 


X fid L 


x gd L 


tan /3 = 30 


tan (3 = 50 


tan/3 = 30 


tan (3 = 50 


1 


1 


0.005 


2.3xl0~ 3 


0.020 


0.009 


0.25 


1 


0.010 


4.7xl0- 3 


0.039 


0.019 


0.25 


4 


0.011 


5.4xl0~ 3 


0.045 


0.022 



Table 2: Present constraints on the off-diagonal squark mass terms which contribute, via 
scalar amplitudes, to AM# d and AMb s ■ All the bounds have been obtained for Ma = 200 
GeV and scale linearly with Ma- 



Notice that the value of the effective coupling \^q[{.5rr{ll))^^ll{rr))^A I nee ded to 
reach the SM level is one order of magnitude smaller than the corresponding one 
required by |Re(5f L )i2|. The precise bounds imposed by Eq. ( |5.5[) on the effective 
coupling |Re[(5^ ij , ( . LL j) 2 3(5f , L ( ij . J? - ) )3i] |, for different values of supersymmetric param- 
eters, are shown in Fig. |^ and reported in Table [p. These bounds indicate that if 
(S^ L )ij and {5ft R )ij are of the same order, we cannot saturate at the same time the 
present limits on B(B S — ► and B{B<i — > On the other hand, this 

possibility remains open if there exists a strong hierarchy between the two types of 
couplings. As we shall see, this conclusion is strongly reinforced by the constraints 
from B-B mixing. 

5.3 B—B mixing 

So far we compared the bounds on the mass insertions relevant to AF = 1 scalar 
amplitudes, at large tan/5, with constraints imposed by AF = 2 vector-type am- 
plitudes (the gluino-box diagrams), which are independent from tan/3. However, as 
pointed out first in |^SJ (see also |29[), at large tan/3 scalar amplitudes could also 
play a significant role in AF = 2 transitions. Indeed contracting the complex 
field in Eq. ( |4.2|) with an appropriate FCNC coupling of H* d [hidden in the Hermi- 
tian conjugate part of Eq. fl4.2p] we obtain AF = 2 effective operators of the type 
d R d 3 L d % L d 3 R whose coupling scales like tan 4 (3. 

If Eij is necessarily proportional to yf, as happens in the chargino case, it is 
easy to realize that these AF = 2 scalar operators are always suppressed by a small 
Yukawa coupling (2/^3) ||28|| . A similar situation occurs in the gluino case, if we allow 
only one type of mass insertion to be different from zero, i.e. {S^ L )sj or (5 R ) R )sj. In 
these cases the only non-trivial constraints are set by B s -B s mixing, where this 
small Yukawa coupling is y s . We have explicitly checked that these bounds are less 
stringent than those set by B(B S — > i.e. they do not exclude the possibility 

that B{B S — > is just below its present experimental limit. 6 

On the other hand, as we learned from the P — > £ + £~ analysis, if we allow 
both LL and RR insertions to be different from zero we can have both E^- and F? 3 



A similar conclusion was recently reported also in Ref. pOl 
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proportional to yt and get rid of any small Yukawa coupling. In this case we obtain 
the following AF = 2 effective Hamiltonian: 



H 



eff (scalar) 
AB=2 



Mlfl+tan^ + ey^)] 2 



'3j 



[l + tan/?(e g + e y y t 2 )] 2 Mj 



b R d 3 L b L d R 



2a a \iM-g 
3vrM 2 

d 



b R d\b L d\ 



R 



X [ x dd L x dd R f( x gd L > x d R d L A)f( x gd R ,Xd L d R A)] > (5-8) 

which have a large potential impact on both AM^ and AMb s ■ Employing the usual 
approximation of degenerate supersymmetric scales, we can write 




200 

W A 

200 



(Sll)3i($RR)31 




\ 50 ) 



[f + 5 (^)l 

/ tan (3 \ 4 
I 50 J 

[§ + 5 (^)] 



4 ! 



4 ' 



(5.9) 
(5.10) 



The bounds obtained by the full analytical expressions, imposing the conditions 
|AM| | < \AM B ^J, are reported in Table @. The comparison between Table |1] and 
Table |2] clearly shows that down-type squark mixing could have a sizeable impact on 
B(B Si d — > /i + /U~) modes only in the presence of a strong hierarchy between LL and 
RR insertions. 

We finally mention that, in principle, one could extract bounds on the combina- 
tion (5^)23(5^)31(5^)23(^^)31 fr° m K-K mixing. We have explicitly checked that 
the bounds imposed by AM K on ^|Re[(5f L ) 2 3(5f L ) 3 i (5^)23(5^)31] I are slightly less 
stringent than those obtained from Kl — > on | Re [(5£ L ) 23(8 RR ) 31] | and reported 

in Table [I]. This happens because, contrary to B s d — > n + n~ modes, we can allow at 
most 0(1) effects in K L — > . 



6. Conclusions 

The rare dilepton decays B S ( i —>■ £ + £~ offer a unique opportunity to study the scalar 
FCNCs which naturally arise in supersymmetry at large tan/3. In this paper we 
have complemented previous studies of these decays modes in supersymmetry 
analysing their sensitivity to non-minimal sources of flavour mixing. 

Similarly to the MFV scenario || |J, we have shown that also in this case it is 
necessary to diagonalize the effective Yukawa interaction — including the d l R d J L Hu 
coupling generated by the new flavour structures — in order to identify the correct 
tan/3 behaviour of scalar FCNCs and to take into account all the leading tan/?- 
enhanced terms. Using this approach, we have corrected/completed previous results 
on A(B S ^ — > £ + £~) in the presence of non-minimal sources of flavour mixing |], [I]]. 
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An important difference with respect to the MFV scenario arises because the 
effective d R d 3 L Hij operator generated by gluino-squark loops is not necessarily pro- 
portional to the Yukawa coupling of the right-handed external quark (yf): allowing 
both left-left and right-right mass insertions to be different from zero, we can gener- 
ate potentially large contributions suppressed only by the bottom Yukawa coupling 
(even for both % ^ 3 and j ^ 3). This leads to extract stringent bounds on the flavour 
off-diagonal down- type squark mixing, at large tan j3, from AM^ a and Kl — > fi + fi~~ . 
These bounds, which have been analysed for the first time as reported in this paper, 
are summarized in Tables 1 and 2. 

As a result of the bounds from AM Bd s and K L — > fi + fi~, we reached the conclu- 
sion that flavour off-diagonal squark mixing can induce order-of-magnitude enhance- 
ments of B(B S) d — > £ + £~) only in the presence of a strong hierarchy between left-left 
and right-right flavour-mixing terms. This scenario, which is apparently fine-tuned, 
could be generated in a natural way in GUT models, as pointed out for instance in 
Ref. [0. If this is the case, one should expect a ratio B(B d — > n + /B(B S — > 
very different from |Kd/Ks| 2 5 i.e. the value expected in the absence of new flavour 
structures. 
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A. Loop functions 

The loop functions appearing in the chargino-squark and gluino-squark penguin 
amplitudes discussed in Sections |3| and can be derived from the recursive for- 
mula [11, 12] 

rl \ I( X i Z li ■ ■ • i z n-2) — I(y, Zl, . . . , Z n -2) / a 1 \ 

I[x,y,z 1 ,...,z n - 2 ) = , (A.l) 



x-y 



starting from the single variable functions 



fix) = — and k(x) = xf(x) . (A. 2) 

x — 1 

Some useful reference values are 

/(M) = ^, /(l, 1, 1) 1 , /(l,l,l,l) = ^. (A.3) 
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